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Aspects of mathematical modeling of the process of melting of lead-containing glasses used for production of 
ceramic paints are discussed. Based on laboratory studies, the stages of physicochemical transformations in 
batch and glass melt are established. The industrial experiments made it possible to determine the optimum 
temperature conditions for production of acid-resistant glasses. This provides the means for justified calcula¬ 
tion of gas consumption in glass melting. 


One of the main criteria of improving the quality of over¬ 
glaze ceramic paints consists in improving their chemical 
stability, primarily, their resistance to the destructive effect of 
acids, as well as preservation of luster in decorated articles, 
which is determined by aesthetic requirements. This can be 
accomplished by improving the quality of the low-melting 
glasses which are the components of such paints. These 
glasses contain lead oxide, which improves various proper¬ 
ties of the product. However, it has a negative effect on the 
acid resistance of glasses and paints made on the basis of 
these glasses. 

Acid resistance depends on glass composition, structure, 
and arrangement of lead cations inside the glass structural 
lattice. 

Research has revealed that it is impossible to completely 
solve the problem of stabilizing the acid resistance of glass 
by varying only the chemical composition of the mixture. A 
determining factor in this case is the temperature/time re¬ 
gime of melting of lead-containing glasses. Such regime en¬ 
sures the preservation of a prescribed lead content and for¬ 
mation of certain structural groups within the glass structure 
whose lead cations exhibit high chemical bond energy. 

Melting of lead-containing glasses is performed in hori¬ 
zontal batch-type rotary furnaces. Crushed and mixed glass 
batch components are charged into the furnace whose tem¬ 
perature is about 1000°C. The process of glass heating and 
melting begins immediately after the charging. 

The contemporary approach to the problem of producing 
low-melting glass is based on mathematical modeling meth¬ 
ods which allow for a sufficiently accurate description of the 
actual processes occurring during batch heating and glass 
melting. 


In the context of the concepts of system analysis based 
on splitting initial problems into subproblems with subse¬ 
quent joining of their solutions, the model of the considered 
problem appears as a combination of the ideal displacement 
model for the gas phase and the ideal mixing model for the 
solid (liquid) phase. The mathematical model equations in 
this case are presented in the following form: 
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where G„ is the gas consumption; C , C s , and C, is the heat 
capacity of the gas, the batch, and the melt; 0 is the gas me- 
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dium temperature in the furnace; L is the furnace length; L, 
<L are the coefficients indicating the existence of convective 
heat transfer to the solid or liquid phases; a p a,, and a, are 
the convective heat transfer coefficients in the gas-solid, gas- 
liquid, and liquid-solid systems; 7j and F s are the tempera¬ 
ture and the surface area of the solid phase; T 2 and F\ are the 
temperature and the surface area of the liquid phase; 7’, and 
F are the temperature and the surface area of the solid SiO, 
particle in the melt; G s and G ] are the quantities of the solid 
and the liquid phases in the system; t is time; P is a coeffi¬ 
cient indicating the existence of the melting process; G™ is 

the quantity of the melting substance in the system; X. is the 
melting heat of the /th component in the system; r a is the 
chemical reaction velocity; AH is the reaction thermal effect; 
Q h is the heat loss to the ambient environment. 

The boundary conditions are T (Z,, r)^ = 0 = 7’ en (x); 

G(i, T)i i=0 = e en (T). 

The initial conditions are T (L, t)| T = 0 = T 0 (L); 

t)| t=o = 0 o (L), 

where T cn and 0 en are the temperatures of the glass melt and 
the gas flow entering the furnace; T (j and 0 O are the tempera¬ 
tures of the glass melt and the gas flow at the initial time mo¬ 
ment. 

Separate calculations for different states of aggregation, 
as well as taking account of heat released in melting of the 
batch components, are carried out by introducing the coeffi¬ 
cients E, 2 , p. which take discrete values [0, 1], 

However, the model in this variant does not take into ac¬ 
count the kinetic specifics of the process and does not spec¬ 
ify which structural and chemical transformations take place 
in the batch and the glass melt under a particular technologi¬ 
cal regime, which is extremely important for comprehending 
the process and constructing a mathematical model which 
would relate the dynamics of the lead cation incorporation 
into the glass structural lattice to the temperature-time condi¬ 
tions for manufacturing acid-resistant glass. Therefore, it 
was necessary to conduct a series of experiments to obtain 
these dependences. The experimental research program con¬ 
templated the following: 

qualitative analysis of the processes which occur in in¬ 
dustrial glass melting; 

X-ray phase analysis of the previous stage results; 

determination of the effect of the batch heating dynamics 
on acid resistance of glass; 

study of the variations in the chemical composition of 
glass melt in the course of industrial glass melting; 

microscope analysis of acid-resistant and acid-nonresis- 
tant crystallized glasses to establish the type of crystal 
syngony; 

X-ray phase analysis of crystallized glasses; 

Glass of the Na 2 0 - B 2 0 - PbO — SiO, systems with 
A1 2 0 3 , ZnO, and TiO, additives was synthesized using in¬ 
dustrial materials (red lead, borax pentahydrate, quartz sand, 
zinc white, alumina, and titanium oxide). The variations in 


the state of aggregation, weight, and color of the batch were 
consecutively analyzed for each 100°C starting with 100°C, 
after 30 min holding at each reference temperature. This ex¬ 
perimental series was conducted up to the temperature of 
1350°C which is the limiting temperature for the production 
technology. The gas medium inside the furnace was kept oxi¬ 
dizing to prevent lead reduction in glass melting. 

At the temperature of 100°C, a certain weight loss was 
registered which continued up to the temperature of 500°C. 
The color of the batch changed within the temperature inter¬ 
val of 300 - 500°C. At the same time, a change in the batch 
volume and apperance of a strong sinter was observed, and at 
the temperature of 700°C the entire surface of the sinter was 
fused while its shape was preserved. An increase in tempera¬ 
ture above 700°C resulted in complete spreading of the melt, 
which formed a smooth surface and still contained the gas¬ 
eous phase and undissolved batch particles, and after holding 
at the final temperature (1350°C) a fluid transparent liquid 
emerged. 

The qualitative changes observed in the course of the ex¬ 
periment were caused by the physical and chemical transfor¬ 
mations occurring inside the initial batch in the glass melting 
process. The weight loss is caused by the removal of chemi¬ 
cally bound water in the dehydration of borax. The change in 
the color of the batch powder is related to the redox pro¬ 
cesses and the chemical reactions occurring in the lead [1], 

The process of phase formation at low heating tempera¬ 
tures was investigated. When the batch is introduced to the 
furnace at a temperature above 100°C, the borax melts and 
decomposes with formation of the compound Na 2 0 • 2B 2 0 3 . 
Assuming that the most fusible eutectic is formed in the 
PbO - B 2 0 3 - Si0 2 system (its melting point is 484°C), on 
further heating (temperature above 100-200°C) 
Na 2 0 • 2B 2 0 5 ought to interact with Pb 3 0 4 and SiO-, and 
form the compounds 5PbO • B,0 3 ■ SiO,, PbO ■ SiO„ and 
5PbO • 4B,0 3 [2], The lead oxide contained in the batch in 
this case should be totally spent on the formation of the 
above compounds. 

An increase in temperature results in the formation of 
compounds in the Na 2 0 - B,0 3 - Si0 2 eutectic system due 
to the existence of Na 2 0 and B,0 3 which are not completely 
converted. The compounds which can be formed in this 
way are Na,0 ■ B,0 3 and Na,0 - 2Si0 2 [2], The remaining 
Na,0 with excessive SiO probably form the compound 
Na,0 • 2SiO, [3], 

The particular scheme of compound formation at a tem¬ 
perature prior to the start of melting is conjectural and takes 
into account the sequence of compound formation with in¬ 
creasing temperature of the eutectics. 

X-ray phase analysis at low batch heating temperatures 
(up to 400°C) only showed the formation of lead silicates 
with different lead oxide contents. The microscope study 
confirmed the presence of these compounds as traces since 
their formation in the solid phase apparently requires pro¬ 
tracted heating at low temperatures. 
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Taking into account all glass-forming components, one 
can expect the formation of more complex eutectic com¬ 
pounds in the Na-,0 - B 2 0 3 - PbO - SiO, system which con¬ 
vert to the vitreous state at the temperature of 500°C with 
complete dissolution of lead, sodium, and boron oxides in 
the glass. The x-ray analysis for this temperature only identi¬ 
fied silica, and the microscope study revealed in addition, 
ZnO, ALOj, and TiO, particles which are not yet dissolved 
in the melt. The x-ray pattern at the temperature of 500°C 
when a consolidated sinter with noticeably altered color is 
formed shows a reflection corresponding to P-SiO-,. This is 
due to the fact that at temperatures close to 500°C the 
eutectic melt in the Na 2 0 - B 2 0 3 - PbO - Si0 2 system is 
formed virtually without any intermediate compounds. 

An increase in temperature causes intense dissolution of 
SiO ? , and zinc, aluminum, and titanium oxides dissolve in 
the melt at temperatures between 900 and 1000°C, which is 
established by the microscope study. A further increase in 
temperature results in decreased viscosity and increased ho¬ 
mogenization of the melt. 

According to the conclusions derived from the experi¬ 
ment and taking into account the high rate of batch heating in 
industrial furnaces and the absence of protracted holding, it 
can be assumed that the formation of a eutectic melt in the 
Na,0 - B 2 0 3 - PbO - Si0 2 system in the glass-melting fur¬ 
nace proceeds virtually without the formation of intermedi¬ 
ate compounds. 

The longevity of a paint coating is ensured by the condi¬ 
tions of the glass melting process. The presence of agitation 
contributes to disappearance of microinhomogeneities in the 
distribution of PbO and B 2 0 3 in the melt. Maintaining the 
maximum glass-melting temperature, which should not lead 
to intense volatilization of PbO, makes it possible to obtain a 
melt without gas phase inclusions and undissolved particles. 
The batch heating rate in industrial furnaces should ensure 
the formation of the strongest bonds of Pb 2+ cations in the 
glass structural lattice, which is accomplished by incorpora¬ 
tion of Pb 2+ cations into vacancies of the glass structural lat¬ 
tice and results in glass compaction [4], 

An experiment in batch heating dynamics was carried 
out in industrial conditions (furnace diameter 1400 mm, fur¬ 
nace length 3200 mm). The experimental data are presented 
in Fig. 1. The main variable parameters were gas consump¬ 
tion, the ratio between gas and air consumption, and the fur¬ 
nace rotational speed. 

A high degree of acid resistance in glass is accomplished 
by keeping the controlling parameters at a certain level to en¬ 
sure the batch heating dynamics in accordance with curve 3 
(Fig. 1). In such conditions the content of the lead cations 
identified in the acetic-acid extract does not exceed 1.7 
mg/dm 2 , which satisfies internationally accepted standards. 
A deviation of the temperature conditions from the optimum 
curve reduces the acid resistance and reduces the luster of 
paints based on these glasses. 



Fig. 1 . Variation of temperature versus the duration of the batch and 
melt heating in glass melting: /, 2 ) accelerated heating rate of batch 
and melt (/ — exceeding the prescribed melting temperature); 3 ) ra¬ 
tional thermal conditions ensuring high acid resistance of glass; 
4 , 5 ) melting temperature below the prescribed temperature) (4 — 
delayed stirring of melt during heating and melting; 5 — absence of 
stirring). 


H PbO’ To 



Fig. 2. Dependence of PbO content in glass on the melting tempera¬ 
ture. 
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Fig. 3. Dynamics of rational gas consumption under the optimum 
glass-melting regime. 

For a deeper understanding of the glass-melting process 
it is expedient to investigate the chemical composition of 
glass melt depending on the technological parameters (fur¬ 
nace length, temperature and duration of melting). 

It is established that the chemical composition of the 
melt varies significantly depending on various factors. This 
is especially true for the lead oxide content which always 
tends to diffuse toward the melt surface and thus decreases 
the surface tension and, accordingly, the surface energy 
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[5, 6], The mass content of lead oxide in samples varies at 
different periods of melting from 30 to 35%, which can pro¬ 
duce increased volatilization. The content of Na 2 0 and B-,0, 
in samples varies to a lesser extent: from 2.5 to 4% and from 
8 to 10%, respectively. 

The variations of the main elements in the glass compo¬ 
sition suggest that increased intensity of glass melt stirring is 
an essential factor of the glass melting process, since it im¬ 
proves the homogeneity of the entire melt body. 

An attempt was made to evaluate the crystallizing prop¬ 
erties of the acid-resistant and acid-nonresistant glass in or¬ 
der to elucidate their structural differences which eventually 
determine the glass quality. Samples of crushed glasses were 
heat-treated for the purpose of their crystallization. 

The microscope analysis of such glasses made it possible 
to determine the extent of their crystallization and the differ¬ 
ence in the crystalline phases they contained. The process of 
crystallization in the acid-resistant glass had been completed, 
as distinct from the acid-nonresistant glass whose degree of 
crystallization was approximately 50%. 

Two phases are distinguished in crystallization of acid- 
resistant glass: one is in cubic syngony (its content is higher) 
and the other in non-cubic (monoclinal) syngony (its content 
is lower). In crystallization of acid-nonresistant glass only 
one phase in the tetragonal modification is revealed. The 
presence of a different crystalline structure point to a possi¬ 
ble difference in the position of lead and boron (which are 
the most important elements determining the prescribed 
glass properties) in the glass structural lattice. 

The x-ray phase analysis data substantiate the different 
variants of crystallization in two types of glass. Lead silicates 
of different modifications are presumably crystallized in the 
acid-resistant glass, whereas lead boron-silicates are crystal¬ 
lized in the acid-nonresistant glass. 


Based on laboratory and industrial experiments, the opti¬ 
mum temperature regime of batch heating and glass melting 
that provides the required acid resistance and luster in the re¬ 
sulting glass was determined (Fig. 1, curve 3). The process 
of incorporation of lead cations in the glass structural lattice 
in this case is implemented in such a way that all parameters 
of the process are efficient at each point of the curve (Fig. 2). 
The natural gas consumption needed to correspond to the 
constructed temperature curve was calculated using the de¬ 
scribed mathematical model (Fig. 3) 

The performed investigations point to the complexity of 
producing lead-containing glass and the strong dependence 
of the process on numerous factors determining the optimum 
properties of glass. The combination of mathematical model¬ 
ing with experimental investigations of all integrated glass¬ 
melting processes makes it possible to determine the opti¬ 
mum production technology parameters to produce glass of 
required quality. 
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